Histone butyrylation in the mouse intestine is mediated by the microbiota and associated with regulation of gene expression.

Post-translational modifications (PTMs) on histones are a key source of regulation on chromatin through impacting cellular processes, including gene expression1. These PTMs often arise from metabolites and are thus impacted by metabolism and environmental cues2-7. One class of metabolically regulated PTMs are histone acylations, which include histone acetylation, butyrylation, crotonylation and propionylation3,8. As these PTMs can be derived from short-chain fatty acids, which are generated by the commensal microbiota in the intestinal lumen9-11, we aimed to define how microbes impact the host intestinal chromatin landscape, mainly in female mice. Here we show that in addition to acetylation, intestinal epithelial cells from the caecum and distal mouse intestine also harbour high levels of butyrylation and propionylation on lysines 9 and 27 of histone H3. We demonstrate that these acylations are regulated by the microbiota and that histone butyrylation is additionally regulated by the metabolite tributyrin. Tributyrin-regulated gene programmes are correlated with histone butyrylation, which is associated with active gene-regulatory elements and levels of gene expression. Together, our study uncovers a regulatory layer of how the microbiota and metabolites influence the intestinal epithelium through chromatin, demonstrating a physiological setting in which histone acylations are dynamically regulated and associated with gene regulation.

Rare subclonal sequencing of breast cancers indicates putative metastatic driver mutations are predominately acquired after dissemination.

BACKGROUND: Evolutionary models of breast cancer progression differ on the extent to which metastatic potential is pre-encoded within primary tumors. Although metastatic recurrences often harbor putative driver mutations that are not detected in their antecedent primary tumor using standard sequencing technologies, whether these mutations were acquired before or after dissemination remains unclear. METHODS: To ascertain whether putative metastatic driver mutations initially deemed specific to the metastasis by whole exome sequencing were, in actuality, present within rare ancestral subclones of the primary tumors from which they arose, we employed error-controlled ultra-deep sequencing (UDS-UMI) coupled with FFPE artifact mitigation by uracil-DNA glycosylase (UDG) to assess the presence of 132 "metastasis-specific" mutations within antecedent primary tumors from 21 patients. Maximum mutation detection sensitivity was ~1% of primary tumor cells. A conceptual framework was developed to estimate relative likelihoods of alternative models of mutation acquisition. RESULTS: The ancestral primary tumor subclone responsible for seeding the metastasis was identified in 29% of patients, implicating several putative drivers in metastatic seeding including LRP5 A65V and PEAK1 K140Q. Despite this, 93% of metastasis-specific mutations in putative metastatic driver genes remained undetected within primary tumors, as did 96% of metastasis-specific mutations in known breast cancer drivers, including ERRB2 V777L, ESR1 D538G, and AKT1 D323H. Strikingly, even in those cases in which the rare ancestral subclone was identified, 87% of metastasis-specific putative driver mutations remained undetected. Modeling indicated that the sequential acquisition of multiple metastasis-specific driver or passenger mutations within the same rare subclonal lineage of the primary tumor was highly improbable. CONCLUSIONS: Our results strongly suggest that metastatic driver mutations are sequentially acquired and selected within the same clonal lineage both before, but more commonly after, dissemination from the primary tumor, and that these mutations are biologically consequential. Despite inherent limitations in sampling archival primary tumors, our findings indicate that tumor cells in most patients continue to undergo clinically relevant genomic evolution after their dissemination from the primary tumor. This provides further evidence that metastatic recurrence is a multi-step, mutation-driven process that extends beyond primary tumor dissemination and underscores the importance of longitudinal tumor assessment to help guide clinical decisions.

Cell-type-specific CAG repeat expansions and toxicity of mutant Huntingtin in human striatum and cerebellum.

Brain region-specific degeneration and somatic expansions of the mutant Huntingtin (mHTT) CAG tract are key features of Huntington's disease (HD). However, the relationships among CAG expansions, death of specific cell types and molecular events associated with these processes are not established. Here, we used fluorescence-activated nuclear sorting (FANS) and deep molecular profiling to gain insight into the properties of cell types of the human striatum and cerebellum in HD and control donors. CAG expansions arise at mHTT in striatal medium spiny neurons (MSNs), cholinergic interneurons and cerebellar Purkinje neurons, and at mutant ATXN3 in MSNs from SCA3 donors. CAG expansions in MSNs are associated with higher levels of MSH2 and MSH3 (forming MutSß), which can inhibit nucleolytic excision of CAG slip-outs by FAN1. Our data support a model in which CAG expansions are necessary but may not be sufficient for cell death and identify transcriptional changes associated with somatic CAG expansions and striatal toxicity.

Selective vulnerability of layer 5a corticostriatal neurons in Huntington's disease.

The properties of the cell types that are selectively vulnerable in Huntington's disease (HD) cortex, the nature of somatic CAG expansions of mHTT in these cells, and their importance in CNS circuitry have not been delineated. Here, we employed serial fluorescence-activated nuclear sorting (sFANS), deep molecular profiling, and single-nucleus RNA sequencing (snRNA-seq) of motor-cortex samples from thirteen predominantly early stage, clinically diagnosed HD donors and selected samples from cingulate, visual, insular, and prefrontal cortices to demonstrate loss of layer 5a pyramidal neurons in HD. Extensive mHTT CAG expansions occur in vulnerable layer 5a pyramidal cells, and in Betz cells, layers 6a and 6b neurons that are resilient in HD. Retrograde tracing experiments in macaque brains identify layer 5a neurons as corticostriatal pyramidal cells. We propose that enhanced somatic mHTT CAG expansion and altered synaptic function act together to cause corticostriatal disconnection and selective neuronal vulnerability in HD cerebral cortex.

B3GALT6 promotes dormant breast cancer cell survival and recurrence by enabling heparan sulfate-mediated FGF signaling.

Breast cancer mortality results from incurable recurrences thought to be seeded by dormant, therapy-refractory residual tumor cells (RTCs). Understanding the mechanisms enabling RTC survival is therefore essential for improving patient outcomes. Here, we derive a dormancy-associated RTC signature that mirrors the transcriptional response to neoadjuvant therapy in patients and is enriched for extracellular matrix-related pathways. In vivo CRISPR-Cas9 screening of dormancy-associated candidate genes identifies the galactosyltransferase B3GALT6 as a functional regulator of RTC fitness. B3GALT6 is required for glycosaminoglycan (GAG) linkage to proteins to generate proteoglycans, and its germline loss of function in patients causes skeletal dysplasias. We find that B3GALT6-mediated biosynthesis of heparan sulfate GAGs predicts poor patient outcomes and promotes tumor recurrence by enhancing dormant RTC survival in multiple contexts, and does so via a B3GALT6-heparan sulfate/HS6ST1-heparan 6-O-sulfation/FGF1-FGFR2 signaling axis. These findings implicate B3GALT6 in cancer and nominate FGFR2 inhibition as a promising approach to eradicate dormant RTCs and prevent recurrence.

Histone bivalency regulates the timing of cerebellar granule cell development.

Developing neurons undergo a progression of morphological and gene expression changes as they transition from neuronal progenitors to mature neurons. Here we used RNA-seq and H3K4me3 and H3K27me3 ChIP-seq to analyze how chromatin modifications control gene expression in a specific type of CNS neuron: the mouse cerebellar granule cell (GC). We found that in proliferating GC progenitors (GCPs), H3K4me3/H3K27me3 bivalency is common at neuronal genes and undergoes dynamic changes that correlate with gene expression during migration and circuit formation. Expressing a fluorescent sensor for bivalent domains revealed subnuclear bivalent foci in proliferating GCPs. Inhibiting H3K27 methyltransferases EZH1 and EZH2 in vitro and in organotypic cerebellar slices dramatically altered the expression of bivalent genes, induced the down-regulation of migration-related genes and up-regulation of synaptic genes, inhibited glial-guided migration, and accelerated terminal differentiation. Thus, histone bivalency is required to regulate the timing of the progression from progenitor cells to mature neurons.

Cell Type Specific CAG Repeat Expansions and Toxicity of Mutant Huntingtin in Human Striatum and Cerebellum.

Brain region-specific degeneration and somatic expansions of the mutant Huntingtin (mHTT) CAG tract are key features of Huntington's disease (HD). However, the relationships between CAG expansions, death of specific cell types, and molecular events associated with these processes are not established. Here we employed fluorescence-activated nuclear sorting (FANS) and deep molecular profiling to gain insight into the properties of cell types of the human striatum and cerebellum in HD and control donors. CAG expansions arise in striatal medium spiny neurons (MSNs) and cholinergic interneurons, in cerebellar Purkinje neurons, and at mATXN3 in MSNs from SCA3 donors. CAG expansions in MSNs are associated with higher levels of MSH2 and MSH3 (forming MutSß), which can inhibit nucleolytic excision of CAG slip-outs by FAN1 in a concentration-dependent manner. Our data indicate that ongoing CAG expansions are not sufficient for cell death, and identify transcriptional changes associated with somatic CAG expansions and striatal toxicity.

Selective Vulnerability of Layer 5a Corticostriatal Neurons in Huntington's Disease.

The properties of the cell types that are selectively vulnerable in Huntington's disease (HD) cortex, the nature of somatic CAG expansions of mHTT in these cells, and their importance in CNS circuitry have not been delineated. Here we employed serial fluorescence activated nuclear sorting (sFANS), deep molecular profiling, and single nucleus RNA sequencing (snRNAseq) to demonstrate that layer 5a pyramidal neurons are vulnerable in primary motor cortex and other cortical areas of HD donors. Extensive mHTT -CAG expansions occur in vulnerable layer 5a pyramidal cells, and in Betz cells, layer 6a, layer 6b neurons that are resilient in HD. Retrograde tracing experiments in macaque brains identify the vulnerable layer 5a neurons as corticostriatal pyramidal cells. We propose that enhanced somatic mHTT -CAG expansion and altered synaptic function act together to cause corticostriatal disconnection and selective neuronal vulnerability in the HD cerebral cortex.

Histone bivalency regulates the timing of cerebellar granule cell development.

Developing neurons undergo a progression of morphological and gene expression changes as they transition from neuronal progenitors to mature, multipolar neurons. Here we use RNA-seq and H3K4me3 and H3K27me3 ChIP-seq to analyze how chromatin modifications control gene expression in a specific type of CNS neuron, the mouse cerebellar granule cell (GC). We find that in proliferating GC progenitors (GCPs), H3K4me3/H3K27me3 bivalency is common at neuronal genes and undergoes dynamic changes that correlate with gene expression during migration and circuit formation. Expressing a fluorescent sensor for bivalent H3K4me3 and H3K27me3 domains revealed subnuclear bivalent foci in proliferating GCPs. Inhibiting H3K27 methyltransferases EZH1 and EZH2 in vitro and in organotypic cerebellar slices dramatically altered the expression of bivalent genes and induced the downregulation of migration-related genes and upregulation of synaptic genes, inhibited glial-guided migration, and accelerated terminal differentiation. Thus, histone bivalency is required to regulate the timing of the progression from progenitor cells to mature neurons.

An allostatic epigenetic memory on chromatin footprints after double-hit acute stress.

Stress induces allostatic responses, whose limits depend on genetic background and the nature of the challenges. Allostatic load reflects the cumulation of these reponses over the course of life. Acute stress is usually associated with adaptive responses, although, depending on the intensity of the stress and individual differences , some may experience maladaptive coping that persists through life and may influence subsequent responses to stressful events, as is the case of post-traumatic stress disorder. We investigated the behavioral traits and epigenetic signatures in a double-hit mouse model of acute stress in which heterotypic stressors (acute swim stress and acute restraint stress) were applied within a 7-day interval period. The ventral hippocampus was isolated to study the footprints of chromatin accessibility driven by exposure to double-hit stress. Using ATAC sequencing to determine regions of open chromatin, we showed that depending on the number of acute stressors, several gene sets related to development, immune function, cell starvation, translation, the cytoskeleton, and DNA modification were reprogrammed in both males and females. Chromatin accessibility for transcription factor binding sites showed that stress altered the accessibility for androgen, glucocorticoid, and mineralocorticoid receptor binding sites (AREs/GREs) at the genome-wide level, with double-hit stressed mice displaying a profile unique from either single hit of acute stress. The investigation of AREs/GREs adjacent to gene coding regions revealed several stress-related genes, including Fkbp5, Zbtb16, and Ddc, whose chromatin accessibility was affected by prior exposure to stress. These data demonstrate that acute stress is not truly acute because it induces allostatic signatures that persist in the epigenome and may manifest when a second challenge hits later in life.

Argonaute-CLIP delineates versatile, functional RNAi networks in Aedes aegypti, a major vector of human viruses.

Argonaute (AGO) proteins bind small RNAs to silence complementary RNA transcripts, and they are central to RNA interference (RNAi). RNAi is critical for regulation of gene expression and antiviral defense in Aedes aegypti mosquitoes, which transmit Zika, chikungunya, dengue, and yellow fever viruses. In mosquitoes, AGO1 mediates miRNA interactions, while AGO2 mediates siRNA interactions. We applied AGO-crosslinking immunoprecipitation (AGO-CLIP) for both AGO1 and AGO2, and we developed a universal software package for CLIP analysis (CLIPflexR), identifying 230 small RNAs and 5,447 small RNA targets that comprise a comprehensive RNAi network map in mosquitoes. RNAi network maps predicted expression levels of small RNA targets in specific tissues. Additionally, this resource identified unexpected, context-dependent AGO2 target preferences, including endogenous viral elements and 3'UTRs. Finally, contrary to current thinking, mosquito AGO2 repressed imperfect targets. These findings expand our understanding of small RNA networks and have broad implications for the study of antiviral RNAi.

High CO2 alters the hypoxia response of the Pacific whiteleg shrimp (Litopenaeus vannamei) transcriptome including known and novel hemocyanin isoforms.

Acclimation to low O2 in many organisms involves changes at the level of the transcriptome. Here we used high-throughput RNA sequencing (RNA-Seq) to explore the global transcriptomic response and specific involvement of a suite of hemocyanin (Hc) subunits to low O2 alone and in combination with high CO2, which naturally co-occurs with low O2. Hepatopancreas mRNA of juvenile L. vannamei exposed to air-saturated water, low O2, or low O2/high CO2 for 4 or 24 h was pooled, sequenced (HiSeq 2500) and assembled (Trinity: 52,190 contigs) to create a deep strand-specific reference transcriptome. Annotation of the assembly revealed sequences encoding the previously described small Hc subunit (HcS), and three full-length isoforms of the large subunit (HcL1-3). In addition to this, a previously unidentified full-length Hc subunit was discovered. Phylogenetic analysis demonstrated the subunit to be a ß-type Hc subunit (denoted HcB), making this the first report of a ß-type hemocyanin subunit in the Penaeoidea. RNAs of individual shrimp were sequenced; regulated genes identified from pairwise comparisons demonstrated a distinct pattern of regulation between prolonged low O2 and low O2/high CO2 treatments by GO term enrichment analysis (Roff-Bentzen, P < 0.0001), showcasing the stabilization of energetically costly translational machinery, mobilization of energy stores, and downregulation of the ubiquitin/proteasomal degradation machinery. Exposure to hypoxia for 24 h resulted in an increase in all of the full-length hemocyanin subunits (HcS, HcL1, HcL2, HcL3, and HcB). The addition of CO2 to hypoxia muted the transcriptomic response of all the Hc subunits to low O2, except for the ß-type subunit.